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1. Introduction 

The quantum filtering equation describes the dynamics of a quantum system undergoing an 
indirect observation continuous in time [1-8]. In the model based on the quantum stochastic 
calculus (QSC) [9-12], worked out by Hudson and Parthasarathy, a role of measuring appa- 
ratus is played by a Bose field [13]. The dynamics of a compound system consisting of the 
quantum system in question and the Bose field is given by a unitary operator. The external 
field is here not only a source of a noise but it also drives the system and carries the infor- 
mation on it [10]. The time development of the posterior state conditioned by a trajectory 
of the results of continuous measurement of the output field is given by the stochastic and 
irreversible equation [5,7,8,14-16]. 

In this paper we apply the theory of measurement continuous in time, based on QSC, to 
derive the quantum filtering equation for a single-mode field in a cavity interacting with the 
Bose field prepared in a coherent state. To get the posterior dynamics of the system we use the 
generating functional approach described in [1,2], but, in contrast to the derivation given by 



Belavkin when the Bose field is initially in the vacuum state, we consider the situation when 
the information about the system is extracted from a measurement of a coherent channel. 
We obtain the nonlinear filtering equation for a direct and heterodyne measurement of the 
signal escaping from cavity. 

The statistics of the output process for the Bose field in a coherent state was considered 
by Barchielli in [17,18], and for a mixture of coherent states by Barchielli and Pero in [19], 
but the estimation of the state of the indirectly observed quantum system was not given 
there. 

2. Quantum stochastic calculus 

In this section we present some basic ideas and notions of QSC in the boson Fock space [9,11], 
which will be needed in our paper. 

We model the electromagnetic field in the symmetric (Bose) Fock space J 7 over the Hilbert 
space /C = C n ® L 2 (IR+) of all quadratically integrable functions from R + into C n , 

jr = C© |jp)/C^ fe j . (1) 

The space /C called 'one-particle space' relates to the states of the photons. For every / G /C 
we define the exponential vector e(/) G T by 

e(f) = (1, /, (2!)- 1 / 2 / ® /, (3!)- 1 / 2 / ® / ® /,...) ■ (2) 

The inner product of exponential vectors in T takes the form 



W)W))t = exp(g\f),c = exp 
The normalized vectors 



w poo 

E/ 

.7 = 1 J ° 



(3) 



*(/) = exp(-^||/|&J e (/), (4) 

are called coherent vectors. In particular the coherent vector t(0) = (1,0,0, . . .) G J 7 is the 
vacuum state. The coherent state t(f) can represent the laser light. 

The linear span T> of all the exponential vectors is dense in J 7 . On the dense domain T> we 
define the annihilation Aj(t), creation Aj(t) and number A^t) processes as follows [9,11]: 

Mt)e(f) = f f^ds e(f) , (5) 
Jo 

(6) 



A](t)e(f) = ± e (f + ex m 



e=0 
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d 

A ij(t)e(f) = -i— e(exp(iAPjjX[o,t))/) 



(7) 



A=0 

where X[o,t) is the indicator function of [0, t), e = (ei, . . . , e n ) G A G K. and (Pijf) k = 8ikfj- 
The operators Aj(t), A^it), Aij(t) satisfy the commutation relations of the form 

[Ai(t),Aj(lf)] = [A\(t),A](t')} = 0, [A tJ (t),A kl (t')] = 5jkAu(t A t') — 8uA kj (t A t') , 

[Mt),A](t')\ = 8^ At', [A 3 (t),A kl (t')} = 5 jk MtAt'), [A kl (t),A](t')\ = 5 l3 A\(tAt'), 

(8) 

where t At' = min(t, t'). 

The Fock space T has a continuous tensor product structure, i.e. 

F = ^*[o,t) ® -Ffcoo) , (9) 

where .T-jo.t) and JF^oo) are the symmetric Fock spaces over C n <g>L 2 ([0, t)) and C n ®L 2 ([t, 00)), 
respectively. Let ft be a Hilbert space associated with some quantum system. The family 
{M(t), t > 0} of operators on ft <E> J 7 is called a quantum adapted process if M(t) acts as 
the identity operator in J-[t,oo) and can act in non-trivial way in ft ® ^o,*) [9]- 

Making use of the factorisation property (9), Hudson and Parthasarathy gave a rigorous 
meaning to the quantum stochastic equations (QSDE) of the type [9, 11] 

n / n \ 

dM(t) = E dA ^ + E iV) dA,-(f) + £>,-(*) dAt(t) + C(f)df , (10) 

i=i \i=i / 

where M(t), Fji(t), Ej(t), Dj(t), C(t) are adapted processes on ft <S> J 7 - The increments 
dA,(t), dA-(t), dyljj(i) acting in ft (g) J-j^t+dt) commute with any adapted process iV(t) in 
ft® J 7 . 

To work with these equations we use the quantum Ito's rule which reads as follows. If 
M'(t) is the process which satisfies an equation of the type (10), then the differential of the 
product M(t)M'(t) is given by the formula [9, 11] 

d(M(t)M'(t)) = dM(t)M'(t) + M(t)dM'(t) + dM(t)dM'(t) , (11) 

where dM(t) dM'(t) can be computed with the help of the multiplication table 

dAi(t)dA](t) = Stjdt, dAi(t)dA kj (t) = S ik dAj(t) , 

dA kj (t)dA\(t) = 5 Jt dAl(t), dA tJ (t)dA kl (t) = 5 jk dA u (t), (12) 
and all other products vanish. 
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3. Filtering equation for photon counting measurement 

Let us consider a single-mode field in a cavity (system S) interacting with the Bose field. 
We assume that the unitary evolution operator U (t) of this compound system satisfies the 
QSDE of the form 

dU(t) = (^bdA\t)-^tfdA(t)-^tfbdt-^_Hd?jU(t), U(0) = I, (13) 

where b stands for an annihilation operator, H = hu (tfb + |) is the hamiltonian of S and 
/i G K is a real coupling constant. The equation (13) is written in the interaction picture with 
respect to the free dynamics of the Bose field [17, 18]. The field described by the operators 
A(t), A'(t), and A(t) is called the input field [10,20] that is the field before the interaction 
with the system S. The output processes, carrying the information about S, are given by 
the formulae: 

A ont (t) = U\t)A(t)U{t) , A out (t) = U\t)A{t)U{t). (14) 
In this section we shall consider a direct detection of the output signal 

A out (t) = ^ (dA(t') + y/jibv cL4t(0 + JjLb\ dA(t') +fibl b t , dt') , (15) 

where b t = W(t)bU(t). The formula (15) can be obtained from (14) and (13) by applying 
the rules of QSC. Note that the self-adjoint process yl out (t) satisfies the non-demolition 
principle [1,2] 

[A out (t'),U\t)ZU(t)} = 0<t'<t (16) 
for any observable Z of the system S, and the commutativity condition 

[A out (t), A ont (t')] = Vf, t'>0. (17) 

Therefore A ont (t) can be represented by a classical random measure on R with values in 
{0,1,2,...}. 

Let us suppose that one realizes a measurement of the output process (15) for the Bose 
field initially prepared in the coherent state t(f)- To get the filtering equation corresponding 
to the observation of a coherent channel, we introduce the generating map, g(k,t), defined 
by 

g(k,t):Z^g(k,t)[Z], 
(i;\g(k,t)[Zty) = ^®LU)\G ou \k,t)Z t ^®LU))i (18) 

where 

t 

G out (k,t) = expjy ln£;(Odyl out (0}, (19) 
o 



4 



Z t = W(t)ZU(t) is the Heisenberg operator of S, ip stands for the initial state of S, and k 
is a complex measurable test function satisfying < \k\ < 1. 

An explicit expression for the generating map (18) can be obtained by solving the differ- 
ential equation for g(k,t). In order to find this equation, one should first obtain the QSDE 
for the operator vr° ut (t,Z) = W(t)G(k,t)ZU(t), where G(k,t) = U(t)G out (k,t)W(t). By the 
rules of QSC, one can easily check that 



dnr(t,Z) = -(Klnr(t,Z) + nr(t,Z)K t )dt + k(t)fib\nr(t,Z)b t dt 

+vWk Ut (*, Z) - (t, Z)b\)dA(t) + y/jHi$*(t, Z)b t - b t n^(t, Z))dA\t) 
+(k(t) - l)(nr\t,Z)dA(t) + ^bUr\t,Z)dA(t) + ^nr%Z)b t dA\t)) , (21) 

where K t = U\t)KU{t) and K — |if + %tfb. 



By taking the mean value of both sides of the equation (21) with respect to the state 
rj — jp (g> <,(/), one gets the differential equation for the mean value of vr^ ut (t, Z): 



(dn%*(t,Z)) = ( V (t)\-(K^ k (t,Z) + 7r k (t,Z)K) + k(t)^7r k (t,Z)b 

+ VJ2(b^ k (t, Z) - 7T fc (t, Zp)f(t) + y^(7T k (t, Z)b - bn k (t, Z))J{t) 

+ (A;(t)-l)(7r fe (t,Z)|/(t)| 2 + ^6t^ (t)Z)/(t) + ^^ (t)Z)6 7^) ) | r?(t))dt) (22 ) 



where ir k (t,Z) = G(k,t)Z and r](t) = U(t)r). Hence, the generating map g(k,t) satisfies the 
differential equation 



dG(k,t) 



(k(t) - l)dA(t)G(k,t). 



(20) 



From (13) and (20), one obtains 



^-g(k,t)[Z] = g(k,t)[- K^Z - ZK - V£(Z&t/(t) +bZf(t)) - Z\f(t)\ 2 



+ k(t)(^Zb + y/jitfZf(t) + y/fiZbf(t) + Z\f(t)\ 2 )] 



(23) 



with the initial condition g(k,0)[Z] = Z. 

The solution to eq. (23) is given by the von Neumann-Dyson series 




(24) 



xS*(t 1 )...S*(t n )Z(t)S(t n )...S(t 1 ), 



(25) 



where 



Z{t) = e - Lt «Ze- L W, 



(26) 




(27) 
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with 

t 

Lit) = Kt + J ^&t /( ^ )+ L/^OTj dt , (28) 



Let r = (ti, t 2 , . . . , t n ) be the trajectory of the observed counting process A ont (t) up to t and 
E* is a set of all finite chains |r| = n e {0, 1,2, . . .}. By introducing the stochastic operator 
V(r 1 1) = e~ L ^S(t n ) . . . S(ti), one can rewrite the series (24) in the form 

g(k,t)[Z] = J k(r)V\r\t)ZV(r\t)dr, (29) 

res* 

n n ^ 

where k(r) = Y\ and dr = Yi The stochastic propagator \/(t)(r) = V(r 1 i) defining 

i=i i=i 

for any trajectory r the posterior state ip(t)(r) = V(r\t)^p of 5, can be represented by the 
Ito chronological multiplicative integral 



= e- L « ^ / dt n / d* n _! . . . / d*i S'(t n ) . . . S'(ti) J] dTV(^) . (30) 

n=o{ { { i=l 

In the last formula 

S'(t) = fMQfyjib + fW-iy-W (31) 

and N(t) is a random variable such that for r G S°° one has dN(t)(r) — 1 if t e r and 
diV(i)(r) = if * £ t. 

By differentiation of (30), we get the linear stochastic equation 

dt/(t) = -(x + V^&V(0 + ^^)^dt+(v^fc + /W-^)^Wd7V(0, 1/(0)=/. 

J32) 

Thus, we draw the conclusion that the posterior unnormalized wave function ip(t) = V(t)ip 
satisfies the linear Belavkin equation of the form 

d$(t) = -( K K + ^f(t)+ l ^^(t)dt+^b + f(t)-I^(t)dN(t), $(p) = i>. 

(33) 

The nonlinear filtering equation preserving the normalization of the posterior wave 
function can be calculated from eq. (33). To derive the differential equation for <p(t) = 
(^(t) |?/>(t)) _1 / 2 V>(t) one has to check that 

d = -(m\^+m)(v^b+f(t))\m)dt 

+ m)\ \(v»tf+wj)(v»b+f(t)) -il \m)dN(t), (34) 
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and insert it into the Taylor expansion of d({ip(t)\ip(t)) 1 ^ 2 ) which yields 

d((m\mr 1/2 ) = 

= (m\mr 1/2 \\ (1*^+2^^)*!®) + & 



+ 



dN(t) , (35) 



"/ \ -1/2 

[p(b% + 2^tRe((b) t f(t)) + \f(t)\ 2 ) -1 
where (.) t = (<fi{t)\{-)<p(t)) . By eqs. (33), (35) and the Ito formula 

dip(t) = d (<^(^)|^(0>- 1/2 ) ^) + Wi)|^(0>- 1/2 d^(t) + d ((^(t)!^))- 1 /^ d ^ (i) (36) 

we obtain the nonlinear filtering equation of the form 

d£(f) = ( - K - ^ fit) + p/2 (b% + v^Re((6) t 7(*)))^(*)dt 

, & + /(*)_ _A ft t )dN(t), ^(0)=^(37) 

Y/x(fotfe) t + 2^Re((fe) t /(t)) + |/(0| 2 / 

If the initial state of the system S is given by a density matrix p, then the posterior 
normalized density matrix pit) satisfies the stochastic equation 

dp(t) = -i[tf,p(t)]dt-^{5t &) p(t)}^ 



( (^jib + f(t))p(t)(VJ^ + f(t)) _ ~ t) 

\p(tfb) t + 2^iRe((b) t f(t)) + \f(t)\ 2 PU 
x (dN(t) - p{b% dt - 2^Re ((b)Jit)) dt - \f(t)\ 2 d?j , (38) 

where {c, d} = cd + dc. The validity of the above equation can be checked by differentiating 

p(t) = mt))m)\. 

By taking the stochastic mean of (38), we obtain the master equation 



d 
dt 

for the prior state 



dt a(t) = -±[H,a(t)]-%{b*bMt)} + iJiba(t)tf + ^jl[b f(t) - f(t) ,a(t)} (39) 



a(t) = (p(t))st, (40) 
since the posterior mean value of dN(t) is given by 

(dN(t))(r) = p(b%dt + 2^Re((b) t W))M+\f(t)\ 2 dt. (41) 

Let us note that for f(t) = the equation (38) takes the form of a quantum filtering 
equation for the reservoir prepared initially in the vacuum state, derived, for instance, in [1]. 



In the model worked out by Hudson and Parthasarathy, the range of frequency of the reservoir 
extends from — oo to +00. A coherent state of the continuum mode distribution, having 
a very narrow spectral width, is here only an analogue of a single-mode laser field and 
f(t) = Aexp(— iut + i0) , A > [21]. In the considered experiment the laser light constitutes 
a coherent signal which stimulates the system S and reaches the detector. We observe the 
interference between the laser light and the light emitted by the system S in the same 
channel. Because of the presence of the terms corresponding to the stimulation of the system 
S, we cannot take the limit |/| — > 00 in (38). 

4. Heterodyne measurement. Transition from counting to diffusion process 

In the heterodyne detection scheme, depicted in Fig. 1, the measured signal is superposed 
with an auxiliary laser field (local oscillator) by means of the beam splitter [17, 18,21]. For 
a lossless beam splitter of transmissivity T, the field reaching detector can be written as 

A mix (t) = VfA™\t)+iVT^TA lo (t). (42) 

The operators Ai (t), A\ Q (t) represent the local oscillator and satisfy the commutation rela- 
tions (8). The auxiliary field does not interact with S and we assume that its initial state is 
given by the coherent vector t(/i )- 

The output generating operator for the ordinary heterodyne measurement is defined as 

t 

G out (M) = (t(/io)|exp{ J ln£;(OdA mix (Oimix(0}^(/io)>, (43) 



where A mix (t) = dA mix (t)/dt. In this way the degrees of freedom of the auxiliary field 
have been eliminated from our description. One can check using (11) and (12) that 
(dA m i x (t)A m i x (t)) 2 = dA mix (t)A mix (t). The assumptions that the field A out (t) is not lost, i.e. 
T — > 1, together with |/i G | — > 00, such that the product (1 — T)|/i Q | 2 := e~ 2 is fixed [10, 16], 
lead to the formula 

t 

G out (k,t) = expjy lnifc(t / )dy out (0} (44) 


with the output counting process Y out (t) of the form 



Y ont (t) = [ dA out (t') + r -^- dA out \t') + ^ dA out (t') + \dt', (45) 

J S EE 



r(t) is a complex function with modulus \r(t)\ = 1. 



8 



Let us note that from the relations (8) it follows that Y(t) = U(t)Y out (t)W(t) commutes 
with Y(t') for any t and t'. Hence, taking into account that [20] U(t)Y out (t' ) = Y{t')U(t) one 
gets 

[Y out (t),Y out (t')] — Vt, t' >0. (46) 

The property (46) allows us to treat the process (45) as a classical one. 

The generating map g(k,t) for the operator G ont (k,t) defined by (44) satisfies the differ- 
ential equation 

t 



^g(M)[£] = g(M) 



- (K + ^f(t) + ^br{t)/e + 1/2 + r{t)/s\ 2 ) Z 

-Z (K + ^f(t) + v^6^)/e + 1/2 |/(0 + r(t)/ £ | 2 ) (47) 
+A;(t) (y^ f + /(0 + K0/ £ ) Z (^6 + /(0 + r(f)/e)' 
with g(k, 0)[Z] = Z. The solution to this equation can be written in the form 

g(k,t)[Z] = J k(K)V\K\t)ZV(K\t)dK, (48) 
where k denotes the counting trajectory of registered photons for heterodyne measurement, 

oo 

fl* — (J {k C [0,0 : \k\ = n}, and 

n=0 

V(n\t) = e- R ®S(t n )...S(t 1 ), (49) 

where 

5(0 = e R W (yjib + fit) + ^ e"^) , (50) 
t 

R(t) = Kt + J ^b^f(t') + ^b^/e + l/2\f(t')+r(t')/e\ 2 yt'. (51) 
o 

Here, the stochastic propagator V(t)(n) = V(k \ t) given by the formula 

t t„ t2 



V(0 = e~ R VJ2 / dt ™ / d*„_i... / dtiS"(t n ) ... S"(ti) JJdy^), 

n=0 i i i, i=l 





where 



S'(t) = e R W (yjib + fit) + ^& - I^j e~ R W , 



(52) 



(53) 



and Y(t) is a random variable such that dY(t)(n) = 1 for t G k and dY(0(/c) = for t k, 
satisfies the equation 

dV(t) = -R(t)V(t)dt+^b + f(t) + r -^--I^)V(t)dY(t), 9(0) = I. (54) 
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Hence, we get for the posterior unnormalized wave function ip(t) = V(t)ip the linear Belavkin 
equation of the form 

d?(t) = -R(t)$(t)dt+ (y]Ib + f(t) + r -^--I^(t)dY(t), ?(0)=^. (55) 

From the equation (55) one can easily derive the nonlinear filtering equation for the normal- 
ized vector (p{t) = (^(t)|^(t))" 1/2 ^(t), 

d<p(t) = -(K + y/»tff(t) + y/»b^/e)ip(t)dt 

+ (/i/2 (b% + ^Re (<&) t (7(0 + K*)/e))) <p(t)dt 

+ ( (/.<&% + 2v ^Re (<6> t (7(0 + 1/(0 + r(t)/e| 2 )" V2 

x (v^& + /(*) + r(0/e) - /) £(0 dF(t) , (56) 



where (.) t = (ip(t)\(.)<p(t)) is the posterior mean value of an operator of the system S. 
Furthermore, the filtering equation for the posterior density matrix corresponding to (56) 
has the form 

dp(t) = (- l -[H,p{t)\ - ^{b%m} + ^mb j + ^[bW)-b j f(t),p(t)})dt 
+ (epbpitp + ^ b (zW) + K0) P(t) + P(t)^ (ef(t) + r(t)) 

-ep{b%p(t) - 2^Re ((6) t (e/(i) + p(0) 
x(edY(t) -e//(6 t 6) t dt-2 v //IRe ((6)* (e/(i) + r(t)) ) dt - e" 1 |e/(*)+r(t) | 2 dt) 
x (e 2 p(b%) t + 2^pZRe (<&) t (e 2 /(t) + e^t))) + \ef(t) + r(t) | 2 ) _1 . (57) 

Of course, if we put r(t) = in the above equations, we get the results obtained in the 
previous section for a direct observation of the output field. Moreover, since the posterior 
mean value of dY(t) is given by the formula 



(dY(f))(«) = (^(b% + 2 y /JIe- 1 Re ((6) t (e/(t) + r(t))) + ef(t)+r(t) 



dt, 
(58) 



the averaging of (57) on the past leads to the master equation of the form (39). 
Now we can consider the linear transformation of the counting process Y(t) [3], 

dW £ (t) = edY(t) - - , (59) 

£ 

leading in the limit e — > to diffusion observation. One can easily check that 

dW £ (t)dW £ (t) = e 2 dY(t) = edW £ (t) + dt, dW £ (t)dt = 0, (60) 
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such that for W(t) = limW £ (t) we obtain the following Ito rules: 

dW(t)dW(t) = dt, dW(t)dt = 0. (61) 
Taking e — >■ one can get from (57) the filtering equation for the diffusion observation 

dp(t) = (- ^[i?,^)] - |{&t fo ,p( t )} + M 6p(0& + + V^[&7(^) -&V(0,p(0])d* 
+ (yjx7(f){b-(b) t )p{t) + ^r(t)p(t) (&t _ (tf) t )) 

x (dW{t) - 2 Re (yji (b) t 7(tj + /(f)f(f)) dt) . (62) 

Let us note that eq. (62) transforms pure states into pure ones. In such a case, for 
p(t) = \lp(t)) {lp(t)\ eq. (62) takes the form 

diffl = [-K + ^bW)-^f{t)+p^) t b-^\{b) t \ 2 )^{t)dt 

+^pW) (b ~ (6)t) m (dW(t) - 2 Re (yji (b) t ^tj + f(t)^ dt) . (63) 

Let us note that if we take f(t) = 0, then all the obtained formulae reduce to the known 
results for the reservoir initially prepared in the vacuum state [3]. Finally it is worth to 
mention that modelling of the evolution of a quantum state conditioned by the results of 
a continuous measurement is the first step towards quantum feedback control [6,8] with a 
coherent source used as a control field [22]. As the choice of the system's operators does 
not interfere with the derivation of the filtering equations for a measurement of a coherent 
channel, the results of the paper can be easily generalized - in the obtained filtering equations, 
the system's hamiltonian and system's coupling operator can be replaced with arbitrary ones. 
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Figure captions 

Fig. 1. The heterodyne setup 
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